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ABSTRACT 


Thermal  radiation  from  refractory  oxide  materials  is 
a  volume,  rather  than  a  surface,  phenomenon.  General  equations 
describing  emissivity,  transmissivity,  and  reflectivity  for 
scattering  and  absorbing  systems  have  been  derived  in  a  form 
that  allows  determination  of  the  necessary  material  constants 
and  calculation  of  emissivity  under  isothermal  and  nan-isothermal 
conditions.  Methods  of  evaluating  conditions  under  which  these 
equations  are  applicable  have  been  derived,  and  a  simple  expres¬ 
sion  developed  for  radiant  energy  transfer  when  conduction  in 
the  solid  is  negligible  o"-  ~  radiant  energy  transfer 

in  the  solid. 

Experimental  ^paracus  has  bean  developed  and  measure¬ 
ments  of  the  isothermal  emissivity  of  several  refractory  oxides 
and  the  transmissivity  of  single  crystals  have  been  made  over, 
the  wavelength  range  1-15  microns  at  temperatures  up  to  1350"C. 
Emissivity  was  found  to  depend  strongly  on  composition,  micro¬ 
structure  and  wavelength.  Different  grades  of  "alumina"  ceramics 
range  from  about  0.2  to  above  0.9  in  the  near  infrared  (<L— ♦ 

Diffuse  transmissivity  measurements  provided  experimental 
values  from  which  scattering  was  calculated  for  each  sample.  Large 
scattering  values  1.  ^00  were  obtained,  however,  emis- 

sivities  calculated  from  these  measurements  agreed  well  with 
experimentally  measured  data.  Results  show  that  the  normal  emis¬ 
sivity  in  these  systems  is  substantially  larger  than  the  hemispher¬ 
ical  emissivity  in  contrast  to  opaque  solids  where  they  are  nearly 
the  same.  \  . 
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1.0 


INTRODUCTION 


Performance  evaluation  of  ceramic,  ceramic-plastic 
or  ceramic-metal  fiber  radomes,  re-entry  bodies  or  high 
temperature  anti-abrasion  coatings  on  a  metal  substrate 
requires  a  precise  knowledge  of  the  thermal  gradients  and 
radiant  thermal  transfer  through  and  from  the  bodies. 

Energy  transfer  through  these  systems  is  by  two  modes: 

(1)  radiation,  (2)  lattice  conduct ivi  try.  There  exist 
complex  interactions  and  perturbations  between  these  trans¬ 
fer  modes  which  cause  large  deviations  from  linear  ther¬ 
mal  gradients  in  areas  near  interfaces.  These  gradients 
may  be  larger  or  smaller  than  would  be  predicted  on  the 
basis  of  lattice  conduction  alone. 

Emissivity  has  been  widely  accepted  as  a  simple 
dimensionless  constant,  and  has  been  used  in  conjunction 
with  the  blackbody  radiation  equation  to  uniquely  express 
the  quantity  and  quality  of  radiant  energy  emitted  by  a 
heated  object.  For  opaque  materials  such  as  metals  or 
graphite  which  possess  large  extinction  coefficients,  the 
emission  process  is  in  essence  a  surface  phenomena,  unaffected 
by  internal  gradients,  and  the  emissivity  concept  for  these 
systems  is  accurate  and  useful.  However,  for  semi-transparent 
materials  such  as  polycrystalline  ceramics  used  for  high 
temperature  applications,  radiant  energy  emission  is  a 
volume  process.  Thermal  gradients  normal  to  the  radiating 
surface  are  almost  always  present,  and  no  obvious  single 
temperature  is  useful  in  characterizing  such  radiating 
systems.  Thus,  the  surface  emissivity  concept  become  ?  fun¬ 
damentally  unsound. 

The  importance  and  significance  of  these  concepts  to 
heat- transfer  am  temperature  measurement  in  translucent 
materials  such  as  refractory  oxides  and  oxide  coatings  is 
not  widely  appreciated.  However,  using  reasonable  optical 
constants  for  dense  sintered  aluminum  oxide,  our  estimates 
indicate  that  differences  as  great  as  one  hundred  percent 
from  isothermal  emissivity  values  may  occur  in  practical 
systems.  For  example,  at  ItO')  C  the  effective  emissivity 
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for  aluminum  oxide  in  still  air  is  already  thirty-five 
percent  higher  than  the  isothermal  value.  This  difference 
becomes  increasingly  important  as  the  temperature  level 
increases  (roughly  proportional  to  T) ,  and  as  the  air 
velocity  increases.  That  is,  these  effects  are  substantial 
ones  of  practical  as  well  as  theoretical  importance. 

McMahon^,  Gordon^,  and  others  ^  ^  have 

developed  the  -'olume  emissivity  concept  for  describing 
the  radiant  transfer  within  and  from  transparent  (i.e. 
non-scattering)  systems  such  as  glass  or  single  crystals. 
However,  most  important  ceramic  and  clastic  systems  con¬ 
tain  either  pores  or  fibers  as  scattering  centers.  As  a 
result,  these  analytical  methods  and  use  of  a  simple  volume 
emissivity  are  unsound.  There  is  not  at  present  any  sound 
general  basis  for  describing  the  thermal  radiation  charac¬ 
teristics  of  semi-transparent  materials  under  non-isothermal 

conditions. 

The  numerical  method  of  Hofctel'  '  and  the  Hamaker 
differential  equations  consider  both  absoiption  and  scattering 
coefficients.  At  present  neither  of  these  systems  has  been 
subjected  to  much  high  temperature  experimental  verification. 
Thus  it  is  necessary  to  (a)  select  an  analytical  system 
sufficiently  powerful  to  solve  complex  and  difficult  problems 
and  which  will  permit  material  constants  to  be  expressed  in 
a  relatively  simple  form,  and  (b)  to  experimentally  verify 
the  analytical  results. 


2.0  THE  PROBLEM 


The  problem  of  heat  transfer  and  radiant  energy 
interaction  with  absorbing  and  scattering  materials  is 
a  complex  one.  Factors  such  as  absorption  and  reemission 
within  the  body  as  well  as  scattering  must  be  taken  into 
account;  in  addition  to  the  linear  thermal  gradients  in 
the  center  of  thick  sections,  boundary  effects  must  be 
considered;  also  there  are  dimensional  effects  and,  less 
frequently,  polarization  and  phase  have  to  be  handled. 

In  general  there  are  twt  methods  of  approach  to  this  problem. 
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Representative  of  one,  the  numerical  approach,  is 
the  Hottel-Cohen  method'®'.  In  this  method  the  sample  being 
considered  is  broken  up  into  small  zones  of  constant  tempera¬ 
ture.  Radiation  exchange  factors  are  determined  for  each 
possible  zone  interchange.  These  radiation  exchange  factors 
are  then  used  in  energy  balance  equations,  one  equation  for 
each  zone  in  the  system.  Heat  transfer  to  the  zone  by  con¬ 
vection,  radiation  or  by  unsteady  state  enthalpy  change  of 
the  zone  are  conside  ed  in  the  energy  balance  equations. 

When  these  are  solved  with  suitable  boundary  conditions, 
the  required  temperature  and  heat  flux  distribution  inside 
the  sample  are  determined.  The  success  of  this  type  of 
calculation  depends  on  the  speed  of  modern  computers.  One 
of  the  advantages  is  that  (by  making  the  zones  small  enough) 
practically  any  desirable  accuracy  can  be  obtained,  pro¬ 
vided  that  sufficient  machine  time  is  available.  A  main 
disadvantage  is  that  a  solution  is  valid  for  only  one 
specific  case;  this  is  useful  when  actual  heat  transfer 
to  some  known  configuration  with  known  optical  properties 
is  desired,  but  is  not  helpful  for  general  design  purposes, 
or  in  understanding  the  mechanisms  involved.  It  does  not 
provide  a  good  basis  for  defining  material  properties  or 
comparing  the  potential  utility  of  different  kinds  of 
materials. 

The  other  general  method  involves  equations  des¬ 
cribing  the  whole  process  in  terms  of  position  and  the 
parameters  at  these  positions.  The  solution  describes 
the  desired  propertv,  such  as  effective  thermal  conduc¬ 
tivity,  in  terms  of  the  necessary  parameters.  The  main 
disadvantage  cf  this  technique  is  that  it  is  extremely 
complicated  to  solve  the  problem  in  the  most  general  case, 
and  the  solution  nay  be  so  complicated  that  a  computer 
is  needed  to  obtain  numerical  values  from  it.  A  practi¬ 
cal  approach  is  to  make  suitable  approximations  which 
simplify  the  general  problem  to  a  manageable  set  of  equa¬ 
tions.  Such  an  approach  was  taken  by  Haniaker^)  who 
assumed  diffuse  radiation  and  divided  the  flux  into  a 
forward  and  backward  component.  Though  the  limitation 
to  diffuse  radiation  simplifies  the  mathematics  a  great 
deal,  it  is  not  a  serious  loss  of  generality  sinm  practi¬ 
cally  all  cases  ot  heat  transfer  of  practical  interest 
involve  mostly  diffuse  radiation. 


Although  other  analytical  methods  are  possible, 
Hamaker ‘ s  equations  appear  to  be  the  most  useful  system 
presently  available  to  characterize  the  thermal  radiation 
behavior  of  both  isothermal  and  non- isothermal  semi¬ 
transparent  solids.  For  non- isothermal  cases,  thermal 
gradients  as  well  as  optical  coefficients  must  be  consid¬ 
ered  in  calculating  radiation  interchange.  This  procedure 
is  not  currently  used  because  of  heat  transfer  specialists' 
greater  familiarity  with  opaque  materials,  and  because  the 
required  optical  data  are  not  generally  available. 

TWO  important  optical  parameters  are  defined  by 
the  Hamaker  relations9:  the  absorption  coefficient  (a) , 
and  the  scattering  coefficient  (s) .  The  former  is  basic¬ 
ally  an  intrinsic  property  that  is  related  to,  but  not 
equal  to  (o')  ,  the  absorption  coefficient  of  single  crystals. 
The  absorption  coefficient  (a)  is  partially  dependent  on 
the  structurally  sensitive  scattering  coefficient. 

In  order  to  be  able  to  calculate  radiant  energy 
emissivity  from  non- isothermal  semi-transparent  materials, 
values  of  (a)  and  (s)  are  needed  as  a  function  of  tempera¬ 
ture  and  wavelength.  The  reflectivity  and  isothermal 
emissivity  of  infinitely  thick  samples  allow  the  determin¬ 
ation  of  the  ratio  (a/s)  but  not  the  individual  values. 

Thus,  present  emissivity  values  cannot  be  converted  into 
the  required  coefficients.  Measurements  of  the  reflecti¬ 
vity,  transmissivity,  or  absorptivity  of  a  given  material 
as  a  function  of  thickness  does  allow  both  constants  to 
be  determined.  In  general,  it  is  most  convenient  and 
accurate  to  measure  transmissivity  for  scattering  systems; 
although  for  very  thin  layers,  completely  diffuse  radiation 
is  required,  and  for  thick  layers  the  accuracy  decreases 
due  to  loss  of  signal  strength  of  the  transmitted  beam. 
Because  radiation  emissivity  differences  between  isother¬ 
mal  and  non-isothermal  specimens  increase  with  increasing 
temperature,  data  for  optical  coefficients  is  most  diffi¬ 
cult  to  obtain  in  the  region  of  greatest  interest.  Since 
(s)  is  nearly  temperature  independent,  this  problem  can 
best  be  overcome  by  elucidating  the  interrelationships 
between  (a)  ,  ( s)  and  (»)  which  are  more  easily  measured 

at  or  extrapolated  to  high  temperature. 
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3.0 


THEORETICAL  INVESTIGATIONS 


The  differential  equations  presented  by  Hamaker 
provide  a  basis  for  developing  analytical  expressions  of  broad 
applicability.  In  contrast  to  some  other  analytical  systems, 
these  derived  expressions  are  general  enough  tr  include  the 
strong  effect  of  scattering  in  addition  to  absorption  effect. 

For  real  systems,  the  relative  simplicity  of  the  equations 
facilitates  a  study  of  mechanisms  and  trends  and  a  determin¬ 
ation  of  design  data.  Although  the  derived  equations  are  not 
entirely  general,  a  number  of  important  cases  may  be  considered 
and  certain  optical  constants  may  be  calculated  for  a  particu¬ 
lar  material.  Other  analytical  systems,  particularly  numerical 
methods,  may  provide  greater  precision  since  some  general 
assumptions  were  made  for  deviation. 

We  have  modified  Hamaker' s  original  differential 
equations  by  introducing  a  correction  for  the  radiant  energy 
density.  This  correction  appears  as  (n2) ,  the  index  of  refrac¬ 
tion  squared,  in  both  the  differential  equations  and  the  gen¬ 
eral  solutions.  An  additional  correction  for  surface  reflection 
is  included  in  the  boundary  conditions. 

Further  development  of  the  Hamaker  relations  to  a 
form  suitable  for  heat  transfer  calculations  depends  on  a 
clear  understanding  of  the  underlying  assumptions  and  the  con¬ 
ditions  under  which  they  can  be  used.  Two  important  assumptions 
involve  diffuse  radiation  impinging  on  the  body  and  the  nature 
of  scattering  produced  by  a  pore  phase.  In  addition,  the 
equations  contain  an  approximation  for  the  conductive  heat  trans¬ 
fer  which  may  impair  accuracy  under  some  conditions.  Methods 
have  been  considered  and  developed  for  evaluation  of  limits 
arising  from  approximations  and  assumptions  inherent  to  this 
system  of  equations.  For  some  of  the  simpler  cases,  new  analy¬ 
tical  relations  have  been  derived. 
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3.1  Further  Development  of  the  Hamaker  Relations 


The  method  used  is  £>aced  on  a  system  originally 
conceived  by  Schuster and  added  to  by  Hamaker (9) ? 
the  notation  used  here  is  essentially  that  of  Hamaker. 

The  basic  method  is  that  of  dividing  the  flux  into  two 
parts:  one  flowing  in  a  positive  direction,  and  the  other 
in  a  negative  direction.  A  set  of  simultaneous  differential 
equations  is  used  to  describe  these  fluxes  and  the  other 
necessary  parameters.  Since  only  a  forward  and  a  backward 
flux  are  considered  this  is  a  one-dimensional  calculation 
and  therefore  has  as  a  basic  assumption  that  the  incident 
radiation  is  diffuse  (i.e.,  the  intensity  is  equal  for  all 
angles  of  incidence) ,  and  that  the  radiation  scattered 
sideways  is  compensated  for  by  an  equal  contribution  from 
neighbouring  parts  of  the  layer  (i.e.,  the  area  investi¬ 
gated  is  either  small  in  cross  section  compared  with  the 
total  illuminated  cross  section  of  the  sample  or  is  large 
compared  to  the  thickness  of  the  sample) .  This  condition 
is  not  a  severe  limitation  since  most  practical  heat  trans¬ 
fer  problems  are  concerned  with  diffuse  radiation. 

The  treatment  for  the  situations  where  temperature 
gradients  are  present  suffers  from  the  further  limitation 
that  only  total  radiation  is  considered  and  therefore  the 
fact  that  the  wavelength  distribution  of  black-body 
emission  changes  with  temperature  is  not  taken  into  account. 
Also  it  is  assumed  that  the  properties  of  the  material 
change  only  gradually.  This  then  implies  the  assumption 
that  the  temperature  gradient  across  the  sample  which  is 
being  measured  is  small.  Practically  all  the  methods  of 
calculation  in  use  today  also  suffer  from  this  limitation 
and  in  practice  there  are  calculation  schemes  which  can  be 
used  to  alleviate  the  problem. 

3.1.1  Isothermal  Layers,  General  Solutions 

The  total  radiant  flux  is  divided  into  two 

parts: 

I  =  The  flux  in  the  direction  of  the  positive  x  axis 

J  =  The  flux  in  the  direction  of  the  negative  x  axis 
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An  absorption  coefficient,  (a) ,  is  defined  by 
requiring  that  (a  I  dx)  be  the  amount  of  the  radiation 
absorbed  from  the  f.lux  (I)  on  passing  through  an  infinitesi¬ 
mal  layer ,  (dx) ;  a  scattering  coefficient  (s) ,  is  similarly 
defined  by  requiring  that  the  flux  scattered  backward 
from  (I)  (and  therefore  added  to  J)  in  an  infinitesimal 
layer  (dx)  is  (s  I  dx)  .  On  passing  through  this  layer 
(I)  will  then  be  diminished  by  the  amount  absorbed  and  the 
amount  scattered,  but  will  be  increased  by  the  flux  lost  by 
scattering  from  (J) ,  or: 


dl/dx  =  -  (a+s)  I  +  sJ 


(1) 


similarly: 


dJ/dx  =■  (a+s)  J  -  si 


(2) 


The  general  solutions  of  these  equations  can  be 
found  by  putting 


I 


+ 


(3) 


J  = 


ox 


e 


+ 


(4) 


only  two  of  the  four  constants  . being  arbitrary. 

The  solutions  (using  the  same  notation  as  Hamaker)  are  then: 

I  =  A  (l-/*0)  e  X  +  B  (1+tfo)  e~  °X  (5) 

J  =  A  (l+/i0)  e  '  X  +  B  (1-ji..)  e“JcX  (6) 


where  _ _ 

^ o  =  n  a  (a+2s) 


(7) 
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=  -J~ a/(a+2s)  =  a0/(a+2s)  (8) 

both  roots  being  taken  with  a  positive  sign.  In  these 
equations  (A)  and  (B)  are  constants  to  be  determined  by 
the  boundary  conditions. 

3.1.2  Isothermal  Layers?  Specific  solutions 
With  No  Diffusing  screen 

One  of  the  cases  for  which  one  would  like 
specific  solutions  is  that  of  a  layer  placed  in  a  beam 
of  diffuse  radiation  where  there  is  reflection  from  both 
internal  and  external  surfaces. 

At  an  interface  where  the  index  of  refraction 
is  increasing,  let  the  reflectivity  equal  (p0).  At  an 
interface  where  the  index  of  refraction  is  decreasing, 
let  the  reflectivity  equal  (p.).  (The  actual  values  of 
these  constants  can  be  calculated  from  the  index  of  refrac¬ 
tion  and  will  be  discussed  in  a  following  section.) 

The  following  nomenclature  will  be  used  (where 
D  is  the  thickness  of  a  layer) : 

I .  =  The  incident  flux  at  x  =  0 

i 

Io  -  The  forward  flux  immediately  below  the  interface 
x  =  0 

Jo  =  The  backward  flux  immediately  below  the  interface 
x  =  0 

JD  ~  The  backward  flux  immediately  below  the  inter¬ 
face  x  -  D 

ID  =  The  forward  flux  immediately  below  the  inter¬ 
face  x  =  D 

There  is  assumed  to  be  no  incident  flux  on  the 
back  surface  x  =  D. 
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Then  the  boundary  conditions  are  that  at  the  front 

surface  (x  =  0) ,  part  (po)  of  the  incident  radiation  (I.) 

is  reflected  back,  and  part  (l-p0)  is  transmitted.  The 

flux  immediately  below  this  interface  (Io)  is  composed 

of  this  flux  (fl-p0]l.)  plus  that  flux  reflected  from 

the  inner  surface  of  x  =  0  or  (p.  Jo). 

r  1 

or  (at  x  =  0)  : 

Io  =  (1-po) I .  +  p .  Jo  (9) 

r  l  i 


At  the  back  surface  (x  =  D)  since  there  is  no  inci¬ 
dent  radiation,  the  only  flux  is  that  part  (p.)  reflected 
from  the  remaining  forward  flux  (1^) . 

or  (at  x  =  D) : 


(10) 


Substituting  in  these  equations  for  Ic,  Jo,  ID,  and 
JD  from  equations  (5)  and  (6) ,  gives: 


A(l-flo)  +  B  (l+f?o  )  =  (1-po)  I±  +  j.±  A(l+/ic)  +  p  B(1-0O) 

(11) 

and 

A(l+^0)  e  +  B(l-^o)  e  ^  =  p^  A(l-^t»)  c  +  ^  B(l+/^)e 

(12) 


These  are  the  equations  to  be  solved  for  the  constants 
(A)  and  (B)  for  these  particular  boundary  conditions.  They 
are  (when  the  exponentials  are  substituted  for  by  hyperbolic 
functions)  : 


9 


I 

I 

\ 

\ 

1^  e  a°D  (1 -pQ)  [Po (1+p^) - (1-pJ  ] 


2  ([/302{I+pJ2  +  (l-p^)2]sinh  a0D  +  2  30  (1-p^2)  cosh  a0D] 

(13) 


It  e°oD  (1-po)  [  0O (1+pJ  +  U-p.)  ] 


2([^o  (1+p.)  +  (1-p.)2]  sinh  CoD  +  20o  (1-p.2)  cosh  OoD] 

3b  X  1 

(14) 


using  these  values  in  equations  (5)  and  (6)  gives  the  follow¬ 
ing  expressions  for  I  (the  forward  flux  at  x)  and  J  (the 
backward  flux  at  x) :  x  A 


I.  <l-po)  t[0o<l+p.)  +  (l-p.)  ](l+|3t,)e'0oXe0oD 

1  r  rl  rx 

+[0o (1+p.) -(1-p.) ] (1-0O) e°oXe~CToD] 

2  ( [0O  2  (1+pJ  2  +  (l-pj^lsinh  a0D  +  20  0  (1-p^2)  cosh  a0  d) 

(15) 


I.  (1-po)  ( [0o (1+p.) +  ( 1-p.) ) (1-Po) e'a°XeCoD 

+  [0O (1+p^) - (1-pJ  ]  (1+0O)  e°oXe  a°D  } 
2 f l P0 2 (1+p . ) 2  +  ( 1-p . ) 2 ] sinh  o0D  +  20o ( 1-p . 2) cosh  0oD) 

r  x  i  i 


(16) 


In  practice,  it  is  impossible  to  check  these  quantities 
experimentally;  what  can  be  checked,  however,  is  the  trans- 
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missivity  and  the  reflectivity.  To  arrive  at  these  qualities 
we  first  determine  the  forward  flux  immediately  under  the 
back  surface  (I  )  by  substituting  (D)  for  (x)  in  equation 
(15). 


Then: 


Ii  2tf  o  (1-po) 


[ tf o 2 ( 1+p . ) 2+ (1-p J 2 ] sinh  OoD  +  2^o (1-p . 2) cosh  0oD 


(17) 


The  transmissivity  (t)  is  then  the  ratio  of  the  amount 

of  radiation  of  the  above  that  gets  through  the  surface 

( [ 1—  p .  ] I_)  to  the  incident  radiation. 
r  l  D 


Or: 


=  I  (l-p.)/l. 
D  r  i  l 


(18) 


giving  for  the  transmissivity: 


2lh  (1-fs,  )  (1-,  .  ) 


[tf,  2(l+..)2+  (l-p.) 2lsinh  o0D  +  2tf0  (l-t.2)cosh  ' ,  D 

r  l  r  l  '  i 


(19) 


The  reflectivity  (,  )  can  be  found  similarly  by  adding  to 
the  fraction  of  incident  radiation  reflected  from  the  front 
surface  (pc.I.)  the  amount  of  backward  flux  that  gets  through 
the  int  erface  (fl-t  .  ]Jc). 


Then: 


r 


(1-,  .)  2-tfo2(l-,  i-2,  ,  )  (1+,  i)  jsinh  '  D+2/'  o  (,  +,  J  (1-,  ^  -osh 


v.2  (1+,  . )  +  ( l-(  . )  2  ’  sinh  '  D+2/V  ( 1-,  .  2)  cosh 


D 


D 


(20) 
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It  is  also  possible  to  calculate  the  absorptivity  (a)  of  the 
layer  since  **  +  p  +  t  =  1.  It  is: 


2&o  (l-p0)  Oo  (1+p.)  sinh  a0D  +  (1-p.)  [cosh  a0D-l]j 
a  =  — - — - — - - - - - 1 - - — - 

f^o2(l+Pi) *  +  (l-pi)2]  sinh  a0D  +  2fi0(l-p.2)  cosh  a0D 


(21) 


This  is  also  the  emission  of  the  layer  relative  to 
black  body  radiation  according  to  Kirchhoff's  law. 


3.1.3  Determining  Optical  Constants  From  Transmissivity 
Measurements 

One  of  the  objects  of  making  transmissivity  measurements 
is  to  use  them  to  calculate  optical  constants  of  the  material. 

In  order  to  do  this,  equation  (19)  for  the  transmissivity  of  the 
material  has  to  be  solved  for  the  constants.  Cross  multiplying 
in  the  equation  (19)  gives: 


2  0O  (l-p0) (l-pi)  = 


2  2  2  2 
t  0O  (1+p^)  sinh  a0D  +  T(l-p^)  sinh  a0D  +  20oT(l-p^  )  cosh  a0D 


(22) 


Or,  by  regrouping  the  terms: 


t  (1-p^)  sinh  o0D  = 


fic  ( 2  [  ( l-p0 )  (1-p . ) -t (1-p . 2) cosh  0oDj -x^o (1+p . ) 2  sinh  o0D]  (23) 

11  x 
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If  we  consider  two  layers  of  thicknesses  (D  )  and 
(D^)and  transmissivities  (i^)  and  (t^)  respectively  then 
(dividing  by  equal  quantities)  : 


t , (1-p . )  sinh  a0D. 

1  r  1  1 


T„(l-n.)  sinh  auD. 

2  r  l  2 


(24) 


2  [  ( 1-po )  (1-pJ  -x^  (1-p^2)  cosh  OoD^  ]-T^o  (1+pJ  2  sinh  o0D 

2 [ ( 1— p o )  (1-pJ -t2 (1-p^2) cosh  OoD^ ]~t2^0 (1+pJ 2  sinh  aDD2 


or  (again  cross  multiplying) : 


3  3 

2  (1-po) (1-pJ  sinh  a oD^- 2i  ^  t2  ( It p  .  )  ( 1— p . )  sinh  cj o D ^  cosh  o  o  2 

2  2 

-P°  t2  (1+pJ  (1-pJ  sinh  OoD^  sinh  0oD^  = 

3  3 

2  t_ (1-po) (1-p . )  sinh  a0D0-2T,  x _ ( 1+p . ) (1-p . )  sinh  o0D  cosh  c0D 
2  r  ri  2  1  2  r  1  r  1  2  1 

2  2 

-PoT.  x _  (1+p.)  (1-p.)  sinh  0oD  sinh  acD 

12^1  r  x  1  2 

(25) 


All  the  terms  involving  (£0)  drop  out  of  the  above 
equation  giving  (having  made  use  of  the  identity  sinh  [x  -y]  e 
sinh  x  cosh  y-cosh  x  sinh  y)  : 

sinh  o0D^  sinh  0oD2  (1+p^)  sinh  oD  (D^-D2) 


T 


2 


(1-  po  ) 


(2c; 
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Or,  if  sample  thicknesses  are  chosen  such  that 


Then: 

sinh  2  '0D 


T 


2 


D  =•  2  D  2D 

-L  ^ 


sinh  o0D 


U+l- ) 

-  sinh  o0D 

(1-po) 


but  sinh  2x  “  2  sinh  x  cosh  x 

and : 


(27) 


(28) 


(29) 


2  sinh  OoD  cosh  o0D  sinh  GoD  (1+p.) 

- _  -  =  — —  sinh  o 0D 

T  T  (1-po)  (30) 

T 2  [  (1+p  J  +  (l-po) ] 

and  cosh  ocD  = - “ 

2  t  (1-po)  (31) 

allowing  one  to  calculate  (a0)  from  two  transmissivity  measure¬ 
ments. 

Once  (g0)  is  known,  (fio)  can  be  found  either  by  solving 
equation  (19)  by  the  quadratic  formula,  or  with  an  electronic 
data  processing  machine. 


3.1.4  Isothermal  Layers  with  Diffusing  Screen 
Corrections. 

In  order  to  approximate  the  theoretical  conditions  in 
the  experimental  measurements,  a  screen  can  be  used  before 
the  sample  to  diffuse  the  radiation.  Since  this  screen  will 
reflect  back  onto  the  sample  some  of  the  radiation  which  the 
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sample  reflects,  a  correccion  must  be  made  for  it.  The  flux 
incident  on  the  sample  is,  however,  the  initial  flux  passing 
through  the  screen  (I  )  enhanced  by  the  radiation  rereflected 
from  the  screen.  If  (I.)  is  the  total  flux  in  this  case  and 
(p  )  is  the  reflectivity  of  the  screen,  (and  the  reflectivity 
of  the  sample  (p)  is  the  same  as  that  calculated  in  the  pre¬ 
vious  section) , 


Then: 


I . 

i 


+ 


2  2 

P  Ps 


(32) 


or 


I . 
l 


2 

n=0 


(33) 


but 

CO 

2  x11  =  l/(l-x)  for  0<x<l  (34) 

n=o 


since  both  0<p<l  and  0<p  <1,  0<pp„<X 

S  o 


and  I.  =  I  /(1-pp  ) 
l  h  s 


(35) 


If  we  introduce  this  result  into  equation  (19)  for  the  trans¬ 
missivity  we  obtain  the  following  result  where  (t)  in  this 
case  is  now  the  ratio  of  the  transmitted  flux  to  the  flux 
measured  from  the  screen  with  the  sample  removed. 


2  fio  (1-po)  (1-pJ 

(1-PP  ) l [tfo2(l+p.) 2+ ;i-p.) 2]sinh  OoD  +  23c (1-p. 2)cosh  ocD] 

S  11  1 

(36) 

If  the  value  of  (p)  from  equation  (20)  is  introduced 
into  equation  (36)  one  obtains: 


2/3  o  (1-po)  (1-p. ) 


I0o2{l+p.)  [l+p.+p  (l-p.-2p0) ]  +(l-p.)2(l-p  )  J  sinh  a0D 

X  1  5  1  X-  S 


*-2/30  (l-p^)  [  1+p^-p^  (po+p^)  ]cosh 


(37) 


This  equation  can  then  be  solved  for  (aQ)  [by  steps 
similar  to  equations  (22)  -  (31) 1: 


sinh  OoD^  sinh  a0  °2  tl+Pi-Ps(po+Pi)  1 

T,  T  (1-po) 


sinh  a0 (D^-D2) 


(38) 


where  (t.)  and  (-r  )  are  the  transmissivities  (as  measured  by 
the  methods  described  here)  of  two  layers  of  thicknesses  (D^) 
and  (D2>  respectively. 

Or  if  thicknesses  are  chosen  such  that: 


D1  =  2D2  5  20 


cosh  o0D  = 


T  2  1  T  1  [  !+Pi-  Pg(po  +  Pi)  HU-po)  } 

2  (1-po) 


(39) 


(40) 


3.1.5  Determination  of  the  Value  of  a0  and  3o  From 

Bxperimen tally  Obtained  Values  of  Transmissivity 


From  equation  (38) ,  0o  can  be  determined  from  measure¬ 
ment  of  the  total  reflectivity  of  a  surface  where  the  index  of 
refraction  is  increasing,  po,  and  the  total  reflectivity  when 
the  index  of  refraction  is  decreasing,  p.,  the  reflectivity 
of  a  diffusing  screen  placed  before  the  sample,  p  ,  and  tne 
transmissivity  and  of  layers  of  thicknesses,  and 
D2  respectively. 
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Since  the  solution  of  this  equation  by  hand  is  laborious 
and  time  consuming,  a  program  has  been  obtained  to  solve  it 
by  machine.  The  Fortran  statements  of  this  program  are 
given  in  Appendix  I.  This  program  was  written  to  be  run  on 
an  IBM  7090  electronic  data  processing  system.  The  other 
material  constant,  the  corrected  emissivity  for  transfer 
between  two  parallel  plates,  Po,  can  then  be  obtained  from 
the  expression  for  t  given  in  equation  (37). 

A  Fortran  program  to  solve  this  equation  for  Po  using 
the  value  of  a0  obtained  above  and  a  value  of  t  and  D 
measured,  and  the  known  values  of  the  other  constants  is 
given  in  Appendix  II.  These  two  programs  allow  one  to  obtain 
the  materials  constants  P0  and  aQ  from  transmissivity  measure¬ 
ments  of  two  or  more  layers  of  different  thicknesses.  These 
programs  have  been  debugged  and  tested  and  the  answers  they 
produce  agree  with  answers  previously  obtained  by  hand. 


3. 1. 6  Surface  Reflections 


In  order  to  use  the  result  derived  here  for  practical 
calculations,  it  is  necessary  to  know  the  value  of  the 
reflectivity  at  the  interfaces.  The  case  of  interest  here 
is  that  of  diffuse  radiation.  Then  the  intensity  is  equal 
for  any  angle  of  incidence,  and  the  flux  is  then  the  intensity 
times  the  cosine  of  the  angle  of  incidence  (as  measured  from 
the  angle  to  the  normal) . 

I r  (e)  is  the  angle  of  incidence,  the  portion  of 
radiation  coming  from  that  part  of  the  hemisphere  between 
(e)  and  (e+  do)  is  sin  ©  do  (times  a  constant  factor  which 
will  be  normalized) .  The  total  flux  at  this  angle  of  inci¬ 
dence  is  then  cos  e  sin  ©  do.  If  p  is  the  reflectivity 
of  the  interface  for  a  ray  at  angle  (e) ,  then  pQ  cos  o 
sin  ©  do  is  the  amount  reflected  for  each  incremental  angle. 

To  find  the  total  reflectivity  then,  one  would  integr^e 
this  expression  and  divide  by  the  total  radiation  (£  sin  © 
cos  e  de) (since  we  are  considering  here  only  an  interface 
where  the  index  of  refraction  is  increasing, this  reflectivity 
is  p0)  : 
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rV/2 

4  p  sin  e  cos  e  de 
Fe 


po  = 


fir/ 2 

•o  sin  0  cos  ©  d  © 


(41) 


The  value  of  (p  )  is  given  by  the  following  Fresnel  equation 
for  reflection  of  randomly  polarized  light  where  (©)  is  the 
angl  '  of  incidence,  (y)  is  the  angle  of  refraction,  (sin  ©  = 
n  sin  y) ,  and  (n)  is  the  index  of  refraction. 


Pe  = 


•  2, 

s in  (e-y) 


sin 


(e+y) 


tan  (e-y) 
tan^ (e+y) 


(42) 


Therefore: 


fir/2  2  2  2  2 

Jo  sin  ©  cos  ©[sin  (©-y)/sin  (©+y)+tan  (e-y) /tan  (e+y)  ]d© 


po  = 


.  rir/2  .  , 

2  /  sm  ©  cos  ©  de 

O 


(43) 


This  integration  has  been  carried  out  by  Walsh 
result  that: 


(12) 


with  the 


po  = 


1 

—  + 
2 


(n-1)  (3n+l) 


6  (n+1) 


[n2(n2-l)2] 
Un2  +  l)3] 


log 


(n-1) 

(n+1) 


2n^ (n^+2n-l) 
(n2+l)  (n4-l) 


[  8n4(n4+l)  ] 

[  (n2+l)  (n4-l>  2  ] 


log  n 


(44) 


18 


which  allows  one  to  calculate  the  reflectivity  of  an  inter¬ 
face  for  diffuse  radiation  as  a  function  of  the  index  of 
refraction.  These  calculations  have  been  carried  out  by 
Ryde  and  Cooper ^^)  with  the  results  shown  in  Table  I?  this 
table  compares  the  total  reflectivity  for  diffuse  radiation 
to  that  at  perpendicular  incidence. 


At  an  interface  where  the  index  of  refraction  is 
decreasing,  the  situation  is  more  complex.  Here  any  radia¬ 
tion  arriving  at  an  angle  to  the  normal  greater  than  sin  "^(l/n) 
is  totally  reflected.  To  find  this  fraction  we  would 
integrate  from  this  angle  to  (tt/2)  and  divide  by  the  total 
flux: 


J  sin 


sin  e  cos  e  de 
-1 ( 1/n) 


total  reflectivity  = 


J  sin  o  cos  9  d© 


(45) 


since 


/  sin  e  cos  e  de  =  ~  sin  e 

2 


(46) 


Total  reflectivity  =  1  - 


n2-l 


n 


(47) 


2  Of  the  flux  that  is  transmitted, (1-reflectivity)  or 

(1/n  ) ,  a  portion  of  this  is  reflected  by  the  Fresnel  reflec¬ 
tion  discussed  above.  This  can  be  shown  to  be  the  same  frac¬ 
tion  as  that  previously  calculated  (p0) .  The  sum  of  these 
two  is  then  the  total  reflectivity  at  a  surface  where  the  index 
of  refraction  is  decreasing  (p^): 


19 


TABLE  1 


Comparison  of  the  total  reflectivity  for 
diffuse  radiation  to  that  at  perpendicular  incidence 
(After  Ryde  and  Cooper (13)) 

iDflfiX  flJL  Reflectivity  for  Reflectivity  for 

Refractio..  Perpendicular  Incidence  Diffuse  Radiation 


1.00 

0 

0 

1.1 

0.0023 

0.026 

1.15 

0.0049 

0.035 

1.2 

0.0083 

0.045 

1.25 

0.012 

0.053 

1.3 

0.017 

0.061 

1.35 

0.022 

0.069 

1.4 

0.028 

0.077 

1.45 

0.033 

0.085 

1.5 

0.040 

0.092 

1.55 

0.047 

0.100 

1.6 

0.053 

0.107 

1.65 

0.060 

0.114 

1.7 

0.067 

0.121 

1.8 

0.082 

0.134 

1.9 

0.096 

0.146 
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+ 


p° 

2 

n 


(48) 


Or: 


n  -1  +  p< 


Pi 


n 


(49) 


This  can  be  an  important  (though  usually  neglected) 
factor  in  heat  transfer  calculations.  For  instance,  for  a 
material  with  an  index  of  refraction  of  1.5,  would  be 
0.595.  For  a  material  of  index  of  refraction  2.0,  p.  would 
be  0.788,  both  factors  being  quite  significant. 

Rearranging  equation  (49) , 


1  —  p  o 


showing  that  the  transmission  through  an  interface  where 
the  index  of  refraction  is  decreasing  is  only  one  divided 
by  n^  of  the  transmission  of  an  interface  where  the  index 
of  refraction  is  increasing. 

3.1.7  Non- I so thermal  Layers 


In  order  to  be  useful  in  heat  transfer  calculations 
this  theory  must  be  extended  to  non-isothermal  situations. 

This  can  be  done  (as  is  shown  by  Hamaker  and  also  by  Schuster) 
by  adding  a  term  in  each  radiation  equation  representing  the 
amount  of  energy  emitted  by  the  infinitesimal  region.  Thij 
is  (e  E  dx)  where  (e)  is  the  emissivity  and  (E)  is  the  black- 
body  radiation  at  the  temperature  at  (x) .  Making  use  of 
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Kirchhoff's  law,  this  term  becomes  (a  E  dx)  where  (a)  is 
the  previously  defined  absorption  coefficient.  An  additional 
heat  balance  equation  is  now  needed  expressing  the  fact 
that  heat  is  neither  accumulated  nor  produced  within  the 
body. 


k  d“T 

~ — —  +  a (I+J)  *  2  a  E  (51) 

dx 


where  (k)  is  the  lattice  thermal  conductivity.  The  first 
term  on  the  left  side  represents  the  heat  accumulated  by 
conduction,*  the  second  term  is  the  heat  absorbed  from  the 
radiation,  and  the  sum  of  these  equals  the  heat  loss  by 
radiation  (the  term  on  the  right) . 

The  total  black-body  radiation  within  a  media  of 
refractive  index  (n)  is  given  by  the  Stefan -Boltzmann  equation 


E 


n 


2  _4 


(52) 


where  (a1)  is  the  Stef an- Boltzmann  radiation  constant  and 
(T)  is  the  absolute  temperature.  If  the  temperature  is 
high  and  the  temperature  gradient  not  too  large  then  (E)  may 
be  represented  by: 


where 


E  =  E0  +  b  (T-T0)  (53) 


b  -  4  a'  n2  Te3  (54) 


and  (T0)  is  a  temperature  close  to  the  actual  temperature, 
and  (E„)  is  the  corresponding  total  radiation.  When  the 
above  equation  holds  the  temperature  may  be  fixed  equally 
as  well  by  (E)  as  by  (T)  and,  since  this  simplifies  matters, 
(fi)  has  been  retained  in  the  equations  rather  than  (T) .  The 
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set  of  simultaneous  differential  equations  is  then: 


dl 

dx 


=  -(a+s)  I  +  sJ  4-  aE 


(55) 


dJ 

dx 

k 

b 


=  (a+s)  J  ~  si-  aE 


d2E 

dx2 


+  a (I+J)  =  2aE 


(56) 


(57) 


Hamaker  shows  that  the  complete  general  solution  of  these 
equations  is: 


I  =  A(l-3)e°X  +  B ( 1+3 )  e  +  C(ax-3)+F 


-ox 


(58) 


J  =  A (1+3) eax  +  B(l-3)e~u~  +  C(ax+3)+F 


-ax 


(59) 


E  = 


.  ax  _  -ax 
-A/ce  -B  k  e  +Cax  +F 


where  a  »  +  2ab  +  a  (a+2s)  =  a0  f  (1+0 

k 


(60) 

(61) 


3  =  a/(a+2s) 


=  2bA(a+2s)  =  2bll/ko 


(62) 

(63) 


And  we  have  introduced  here  the  proper  (n  )  term  which  does 
not  appear  in  Hamaker 's  work. 

To  illustrate  how  this  theory  might  be  used,  we  will 
derive  the  particular  solutions  for  a  layer  receiving  radia¬ 
tion  at  both  surfaces,  and  where  heat  is  being  conducted  away 
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from  the  surfaces.  The  amount  of  heat  being  conducted  away 
from  the  surface  must  equal  that  conducted  to  the  surface 
in  the  solid  giving  one  boundary  condition  at  each  surface. 
The  other  two  boundary  conditions  are  supplied  by  the  radia¬ 
tion  interchange  at  the  surface.  The  temperature  (particu¬ 
larly  at  the  surfaces)  and  the  emitted  fluxes  will  be  solved 
for. 

Using  the  same  notation  as  the  previous  sections , 
immediately  below  the  front  surface  x=0 ,  the  forward  flux 
(Io)  is  equal  to  that  part  of  the  incident  flux  (I.)  which 
isn't  reflected  ([l-p0]l^)  plus  the  amount  of  the  1 

backward  flux  at  this  surface  (J0)  which  was  reflected 

( P  O  )  • 

Therefore; 

Io  «  (1-po)  I.  +  p±Jo  (64) 


Or  substituting  from  equations  (58)  and  (59) : 

A  (1-3)  +  B(l+3)  -CJ3  +  F  =  (l-p0)Ii  +  p  A(l+3)  +  p  B(  1-/3) 

+  P±  c3  +  FPi 

(65) 


Similarly,  immediately  below  the  back  surface  x=D, 
the  backward  flux  (JQ)  is  composed  of  the  part  of  the  inci¬ 
dent  flux  on  this  surface  (J.)  which  is  trfmsmitted  (fl-p0]J.) 

plus  the  part  of  the  forward  flux  at  this  surface  which  is 
reflected  (p^Ip) : 

And: 

JD  =  (l-po)Ji  +  p  ID  (66) 
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or: 

A{l+/3)e°D  +  B ( 1-4*)  e~  D  +  C ( aD+fi)  +F  = 


OD 


(l-p0)Ji+pi  A(l-^)e  +B  H 


(1+^) 


-oD 


P. 


C(aD-&) + 


P 


.F 

1 


(67) 


If  we  define  (t))  as  being  equal  to  the  gradient  at  the  sur¬ 
face  times  (b) : 


dT 

b(3? 


surface 


dE 

dx^  surface 


(68) 


Or,  if  the  heat  is  conducted  away  by  a  gas: 


dE  -b  0 

■n  =  ( - )  =  - - 2. 

dx  surface  k 


(69) 


where  (Qg)  is  the  heat  being  conducted  (or  convected)  away 
by  the  gas. 


Then  since  (by  differentiating  equation  (60)  )  : 


dE 

dx 


*  ox  ,  „  -ax  ,  _ 

A  k  a  e  +  B  a.  a  e  +  C  a  (70) 


The  other  two  boundary  conditions  are: 


1 


dx 


o 


T)  =  -A  /ca-tBA.a  +  Ca 


(71) 

(72) 


and 


(^) 

ax' 


(73) 
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-A 


OD  JL 

a  e  + 


B  k.  a 


e"aD  +  C 


a 


(74) 


The  four  simultaneous  equations  (65)  ,  (67)  ,  (72) 
and  (74)  are  then  solved  for  the  constants  A,  B,  C,  and  P 
where  it  has  been  found  convenient  to  define  a  function  con¬ 
sisting  of  the  denominator: 

let: 

2{l~p  J  (cosh  cD-i)  +  [  ( 1+p^)  (l+/c)+/c  crD(l-p  )  jsinh  oDsetc  (75) 

then: 

(e~aD-l)  {a  (l-p0)  (I±-J^)  +  r)[20(l+p  )+aD(l-p  )  ]} 

A  =  - - - - — 

2  a  etc  (76) 


(e°D-l)  {a (l-p0)  (I.-Ji)+T1[2/3(l+pi)+aD(l-pi)  ]) 

B  ~  1 . . . ’  1 1  "  ~ 

2  c  etc  (77) 

-2  /co  (l-p0)  (I^~  Th  )  sinh  oD+4t)  [  (1-P^)  (cosh  oD-1) 

+  P  (1+p . )  sinh  aD  ] 

c  = - _i__ - 

2  a  etc  (78) 

2on2Ii(  (1-Pi)  (cosh  a  D-l)  +  [  (/3  (1+Pi)  (l+zcJ+zcaDd-pJ  ] sinh  aD] 

+  2an2  J.(  (1-Pi)  (cosh  aD-l)+0(l+p  )  (l+/c)sinh  aD] 

-2arjD(  (1-p  . )  (coshoD-X)  (1+p .)  sinh  aD) 

F  = - - - - - - - 

2  a  etc 

(79) 
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By  introducing  these  constants  into  equations  (58) , 
(59) ,  and  (60)  it  is  now  possible  to  find  the  fluxes  and 
temperature  at  any  point  in  terms  of  the  radiant  and  thermal 
fluxes. 

They  are: 


I 


a(l-p0)  (I.-J.)  (eax(e  °D-1)  (l-0)+(e  -1)  (l+0)e  °x-2k,  (ax-3)  sinh  aD) 

1  l  l 

+ri  ( f  20 (1+p . ) +aD(l~p . )  ][  e°X  (e-oD-l)(  1-0) +e-0X  (e°D-l)  (1+0)  ] 

+4  (ax-0)  [0(l+p^)  sinh  aD  +  (1-p^)  (cosh  aD-1)  ]} 


2  a  etc 


+F  (SO) 


a 


(1-po) (I.-J.)  (eax(e_aD-l) (l+0)+e"ax(eaD-l) (1-0) -2* (ax+0) sinh  aD] 


l  l 


+il  ( [  20  ( l+p±)  +aD ( l-p±)  ][eax(e~aD-l)  (1+0) +e  °x(eUiJ-l)  (1-0)  ] 


-ax ,  aD 


J  = 


+4  (ax+0)  [0  (1+pJ  sinh  aD  +  (1-pJ  (cosh  aD-i]} 


2  a  etc 


+F 


(81) 


■k.  a  (1-po)  (I^-Jj  (eax  (e  °D-l)+e  ux  (e^-1 1  +2ox  sinh  aDj 


-ax,  aD 


■T)U[20(l+p.)  +  aD(l-p.)  ]leax(e"aD-l)+e~ax(eaD-l)  ] 


-4ax  [ 0 ( 1+p^)  sinh  aD  +  ( 1 - P ^ )  (cosh  aD-1) ] j 


2  a  etc 


+F 


(82) 
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These  are  again  not  measurable  quantities.  The  quan¬ 
tities  desired  are  the  fluxes  emitted  at  each  surface  and 
the  temperatures  at  the  surfaces.  At  the  back  of  the  layer, 
(x^D) ,  the  flux  emitted  in  the  forward  direction  (here 
denoted  by  Ie)  is  equal  to  the  fraction  of  the  forward  flux 
immediately  under  this  surface  (ID)  which  is  not  reflected 
at  this  surface  ([l-pillD^  plus  tne  fraction  of  the  incident 
radiation  on  this  surface  which  is  reflected  into  the  for¬ 
ward  direction  (poJ.). 

r  r 


I  =  (l-p.)I_  t  p0J. 
e  ri  D  r  i 


(83) 


Similarly  (where  J  is  the  flux  emitted  in  the  back¬ 
ward  direction  at  the  front  surface) : 


J  *  (1-p.)  Jo  +  pol.  (84) 

e  i  'i 

Finally  the  energy  equivalent  of  the  temperature  at 
the  surfaces  is  found  by  substituting  x=0  and  x=D  into 
equation  (82) .  Actually  in  the  constant  (b),  the  (n^)  term 
should  be  the  index  of  refraction  of  the  material  in  which 
the  particular  quantity  is  measured.  In  order  to  keep  the 
notation  consistent  in  this  part,  the  (n2)  term  will  be  kept 
in  the  constant  (b)  but  the  energy  equivalent  of  temperature 
measured  outside  the  sample  (here  denoted  by  E  and  E  D)  will 
be  divided  by  (n2)  so  that  the  numerical  results  will  Be 
correct.  Before  setting  down  the  results  it  is  desirable 
to  define  the  following  functions  since  most  of  the  equations 
are  symmetrical. 

Let  etc.  =  2(l-p^)  (cosh  oD-1)  +  [  2$  (1+p  J  (l+/c)  +/caD(l-p  J  ]sinh  oD 

(85) 
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2/3(l-p0)  (l+/c)  sinh  oD 


(86) 


2(l-pJ  (cosh  aD-1)  +  [  2/3  (p0+p^)  (1+k.)  +/caD(l-p^)  jsinh  aD 


f2  " 


(87) 


(1-k)  (1-p^)  (cosh  aD-1) + [0 (1+p^)  (l+/c) +/caD(l-p^)  ]sinh 


(88) 


(l+/c)  [  (1-p . )  (cosh  aD-1)  +  /?(l+p.)  sinh  aD] 

i  i 


f4" 


(89) 


(1-pJ  [-2 (cosh  aD-1)  +  aD  sinh  aD] 


f  = 
5 


(90) 


[  aD  (l+/c)  (1-p  .) +2^<c  (1+p  .)]  (cosh  aD-l)+  PoD  ( 1-p  .)  sinh  oD 


n2  etc 


(91) 


It  is  useful  to  note  that, 


fl  +  f2  =  f3  +  f4  =  1 


(92) 
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f4  "  ‘  f5 


(93) 


f 


3 


Then  the  desired  terms  are  (where  £g  is  the 
heat  removed  by  conduction,  and  therefore  t) 


negative  of  the 

=  b  Q  A)  : 

-g 


i 

e 


f,  I.  +  f,  J.  +  f,  <  0 
1  l  2  i  5  -g 


(94) 


J 

e 


E 

ao 


f0  I.  +  f.  J.  -  f_  k  Q 
2  i  1  i  5  -g 


f,  I.  +  f.  J.  -  f_  Q 

3  i  4  a.  6  2 p  -g 


K. 


f  I  +  f.  J.  +  f_  « 

4  i  3  l  6  2&  -9 


(95) 

(96) 

(97) 


This  gives  an  idea  of  how  the  theory  can  be  used  to 
calculate  radiation  transfer;  the  actual  choice  of  parameters 
is  arbitrary.  For  instance,  it  would  be  possible  to  specify 
the  surface  temperatures  and  perhaps  assume  no  incident  fluxes 
(they  might  conveniently  be  chosen  as  being  zero) ,  then  solve 
for  the  emitted  fluxes  and  the  surface  gradients  necessary  to 
maintain  the  given  situation. 


3.1.8  Derivation  of  the  Effective  Thermal  Conductivity 
of  a  Powder. 


Using  the  relations  derived  above,  it  is  possible  to 
calculate  an  effective  thermal  conductivity  for  a  powder  if 
it  is  assumed  that  the  powder  consists  of  layers  of  material 
the  thickness  of  which  is  the  same  as  the  average  particle 
size  of  the  particles  of  the  powder  and  if  the  porosity  of 
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the  model  is  the  same  as  the  porosity  of  the  powder.  An 
amount  of  radiation  is  incident  on  the  first  layer  of 

the  system;  the  radiation  being  emitted  then  from  its 
second  surface  is  Ie^„  This  latter  amount  is  now  incident 
on  the  second  layer  and  is  therefore  the  same  as  I^.  On 
the  last  layer  of  the  system,  the  j-tll  layer,  the  amount  of 
radiation  incident  is  (it  is  in  the  backward  direction) , 

and  this  layer  emits  from  its  other  surface  an  amount  of 
radiation  Je j ,  which  is  then  incident  on  the  next  to  last 
particle  and  would  then  have  the  notation  there  of  J, ^ . 

In  this  manner,  all  the  radiation  terms  can  be  found'  3 
as  functions  of  the  radiation  incident  on  the  first  and  last 
layers.  It  is  assumed  that  the  first  layer  is  the  hottest 
one  and  heat  is  transferred  down  the  stack  to  the  j.^h  layer. 

Since  the  emitted  radiation  from  one  layer  is  inci¬ 
dent  on  the  next,  one  can  transform  equations  (94)  and  (95)  to 
the  following; 


I  -  f.  I  .  -  f-  <  0  *  0 

en  1  e(n-l)  2  e(n+l)  5  -g 


(98) 


and: 


;  .  -  f.  J  ,  ...  +  fc  /ck  Q  =  0 

e(n-l)  1  e(n+l)  5  g  -g 


(99) 


But  since  f^  =  l-f2,  substituting  -Ie  +f2Ie  (n-1) 

for  •flIel«-l|  in  e<?ua*'ion  (98)  and  -Je  (n+1)  +  Ve(n+’ )  for 

-f,J  /  ,  \  in  equation  (99)  ,  the  following  is  obtained: 

1  evn-lj 


I  -I  ,  ,v+  f.[l  ,  .,-J  ,  , .  x  ]~f  -  ^  Q  =  0 

en  e(n-l)  2  e(n-l)  e(n+l)  -g 


(100) 


I 


l 


and: 


J  -J 
en 


e(n+l) 


^^etn-l)  "Je(n+1)  ^ 


0 


Similarly,  substituting  for  £2  ^e(n-M)  in  equation 
^2Ie(n  1)  *n  e<3uat^on  (99)  .  The  following  is  obtained: 


'en  ~  Je(n+1)  ^l^efn-l) 


Je(n+1)  -1  +f! 


Q  = 
-g 


and: 


J 

en 


Ie(n-1)  +  ^l^etn-l) 


Je (n+1) 


-f. 


/c 


Q 

-g 


0 


Applying  equation  (95)  to  the  (n+-J  layer: 


1  £l  ^5 

en  f2  Je(n+1)  f2  Je(n+2)  +  f2  K  -g 


and  equation  (94)  to  the  (n-1)  layer: 


j  s  *  i  (  ^  -  ”r  J  /  m 

en  f2  e(n-l)  f2  e(n-2) 


Q 

f2  -g 


substituting  equation  (7.04)  into  equation  (102)  : 


(101) 


(98)  , 


.02) 


(103) 


(104) 


(105) 
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^e(n+l)  ^l^e(n+2)  +  ^5  "  -g  ^2  Je(n+1) 

tl  f 2  ^  Ie  (n-1)  ~  Je  (n+1)  ^~f2  r5  -g  =  0  (106) 

Or,  since  l~f2  ~  f ^ ,  and  dividing  by  f^: 

J  (  xn  "  J  /■  “  f-,  f*  ,  -  J  ,  ]  +  fc  /c  Q  =  0  (107) 

e(n+l)  e(n+2)  2  e(n-l)  e(n+l)  5  -g 

Similarly,  substituting  equation  (105)  into  equation  (103): 
Ie(n-1)  *  Ie(n-2)  +  f2®-Ie(n~l)  "  Je(n+1)  ■*  f5  <  Q  =0  (108) 

-g 

From  equations  (100) ,  (101) ,  (107) ,  and  (108) ,  it  can 

be  seen  that: 

Jen  ^e(n+l)  Je(n+1)  Je(n--2)  ^e(n-2)  ^(n-l)  = 

^e(n-l)  ^en  ^en  *ei  n+lA  ^e(n+l)  * e(n+2 )  (109) 

Furthermore , 


33 


Xe(:i-l)  "  Je(n+1)  “  Ie(n-2)  ‘  f2  f  Xe  (n-1)  ~Je  (n+1)  ^+f5  *  -g 


Je(n+1)  ^2  ^©(n-l)  ~e (n+1)  J  ’  *5 


-g 


I  /  “  J  /  , '  2  fjl  .  -j  .  .  ] 

e(n+2)  e  (n+1  2  a  (n+1)  e(n+l) 


+  2  f  «;  Q 
b  ~g 


(110) 


Or  by  repeated  substitutions  similar  to  aquation  (109)  ,  if 
there  are  j  particles. 


Ie(n-1)  Je(n+1) 


13  *  Jij  -«-«  f-  [I 


2  ,'',‘e(n-l)  Je)n+1)^ 


+  ( j  —  1)  f.  /c  Q 
5  -g 


(111) 


and: 


Ie(n-1)  Je(n+1) 


I..  -  J. .  +  (j-1)  £  k  Q 
xl _ li _ 5  -  q 

1  +  f2  (j-1) 


(112) 


and: 


Jen  Je(n+1)  Je(n+1)  Je(n+2)  ^en  Je(n+1) 


Ie(n+1)  Ie(n+2) 


f0[l..  -  J.  .  ]-fc  *  Q 
2  ll  1}  5  -g 

1  +  f2  (j-1) 


013) 
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Since  the  amount  of  heat  transferred  by  gas  con¬ 
duction  is  the  same  everywhere  in  the  system,  it  would 
simplify  matters  if  the  following  substitution  is  made. 


(dT)  _  ^Ead(n-1)  Eaon^ 

(dx) g  ”  b  D 

1  P 


The  notation  b^  (which  is  defined  as  being  equal 
to  4  a '  T^)  is  here  used  to  denote  transfer  in  air  where  the 
index  of  refraction  is  unity  as  differentiated  from  b  (which 
equals  4  a'n  T  )  which  is  used  for  transfer  within  a  layer 
of  index  of  refraction  n. 

Where  it  should  be  remembered  that: 


k 

g  (<2x)g 


(  llc  i 


D  represents  the  distance  across  the  pores  and 

E  -E 

ad(n-l)  aon  is  the  temperature  drop  across  the  pore.  The 


latter  can  be  evaluated  from  equations  (96) and  (97) : 


(E 


ad (n-1) 


-E  ] 

aon 


Ie(n-2) 

3  Ie(n-1) 


J  + 
en 


^4  Je(n+1) 


ld£L 

(dx)  g 


e  ( n-  2) 


end  J 


(116) 


Or,  on  substituting  for  I 


en 


from  equations  (lUB)and 


(101) ,  and  using  the  relationship  f^-f^  -  k 


[E  ,,  1N-E  1  = 
ad(n-l)  aon 


[le(n^l)-Je(n+l)][f2-'  f5] 

k  fC 

1  +  K  -3-  [  _  f  ] 

b.D  1  0  x5J 
1  p 


(117) 


Finally,  by  substituting  for  x)“Je(n+l)  ^rom 

equation  (112),  again  solving  for  Ead^n  x) ~Eaon  and  reducing 
to  simplest  terms  we  finally  obtain  the  latter  as  a  function 
of  I . ,  J.  .  and  k  : 

il  13  g 


[E 


ad (n-1) ~  aon 


]  - 


.  "n-yy 


l+-f2  ( j  *,  b  D 


1  P 


5fc..f5(j-l)+l]) 
(118) 


Introducing  equation  (1.18)  into  equation  (113),  the  following 
is  obtained: 

Jen  Je(n+1)  Je(n+1)  Je(n+2)  *on  ^(n+l)  Ie(n+1)  Ie(n+2) 


/ck  f0f,  o 

fcu-V  (v  b^9  -  *  f5  ” 

l  p 

/ck  £,. 

l+f2(j-l)+  b—g  {  -j  [l+f2(j-l)  ]-f5[/cf5(j-l)+l]) 


1  P 


(119) 


Using  the  above  equation,  any  desired  term  can  now  be 
derived.  For  instance: 
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1  1=  1  ! 
el  xl 


—  IP.  —  .  - 

/ck  t. 


l+f2(j-l)+  h~jf  f-jcf  [l+f2(j-l)-f5[/cf5(j-l)+l])  (120) 

Ip 


and: 


n[1irJijHf2 


I  «  X.. 
en  ll 


[  f2f6  -  <c  f  _2]) 

Vo  <» 


similarly: 


^clc  f  "f” 

(j-n+1)  [XjL1-Ji  ]  {f2  +  *  f5  ^ 

1  p 


J  =J. . + 
en  13 


/ck 


/ck  f 

l+f2(j-l)+  b-^2  -j  Cl+f2(j-l)-f  [/cf  (j-l)  +  l]} 

ip 


(121) 


l+f2(j-l)+  h~*  (  £§  [l+f2(j-l)-f5[/c  f5(j-l)+l]) 

1  p  /3 


(122) 


And  substituting  into  equations  (96)  and  (97) : 


2  =  f,I. ,  +  f .J. 

aon  3  xl  4  13  zck  f  f  f 

[lil-J  b~fl^  — fs2]}+ 

1  p 


Kk  f 

l+f2(j-l)+b-~  {  ~f  [l+f2(j-l)  ]-f5(j-l)+l]) 

1  p 


(123) 


and: 
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E  .  =  f„I.,  +  . 

adn  4  ll  3  lj 


£  £  £ 

riil+Ji jH[f3(j-l)-n+lHf2  +  b-^  [  -P  -«fs2])^f£2-«f5]} 

1  P 


/ck 


1  +  f  (j-l)  +hp  (  [l+f0(j-l)  ]-fg[/c  f  ( j-l)  +1] } 

ip 


(124) 


These  relationships  allow  the  calculation  of  any  desired 
radiation  tern  or  temperature  in  a  system  if  the  incident  radia¬ 
tion,  the  conductivity  of  the  gas,  the  scattering  and  absorption 
coef ficients ,  the  average  temperature,  the  conductivity  of  the 
solid,  and  the  average  size  and  number  of  particles  in  the  sys¬ 
tem  are  known. 

An  effective  thermal  conductivity  ke  is  defined  by 
the  following  equation: 


or: 


(Ax) 


(125) 


(126) 


where  Q  is  the  heat  flowing  across  a  unit  area  under  the 
influence  of  a  temperature  drop,  AT,  which  occurs  over  a 
representative  distance,  l\x.  The  length  from  the  top  sur¬ 
face  of  one  layer  to  the  fop  surface  of  the  next  layer  will 
be  taken  as  a  representative  distance.  Then: 


Ax  =  (D  +  D  )  (127) 

P 

where  D  is  the  thickness  of  a  layer,  and  is  the  thickness 
of  a  pore.  The  temperature  drop  across  this  section  is  then: 
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\E  ,  ,  ,.-E  ] 

ad(n-l)  aon 


AT  = 


(128) 


which  from  equation  (123)  is  found  to  be: 


AT  = 


-  3-1  i.l 


/ck_  f  ~±,,  „ 

r  _  .  9  ,26  _  2.  , 

1  f2  +b,D  (_ jr  -  K  f5  >  1 

_ _L_P_ 


f  1+f  (j-l)+—S  (-1  [1+f  (j-1)  ]-f  Uf  (j-D+j) 
i  2  b  D  ^  2  5  5  j 

1  P 


(129) 


The  total  heat  flow  can  be  found  if  any  surface  is 
considered.  There  the  heat  flow  is  equal  to  the  difference 
between  the  forward  and  backward  radiation  fluxes  plus  the 
heat  conducted  by  the  gas.  The  latter  is  equal  to  the  tempera¬ 
ture  drop  across  a  pore  times  the  thermal  conductivity  of  the 
gas  divided  by  the  distance  across  the  pore  or: 


Q  =  I .  -  J  +  r-^r-  [E  .  -E  ] 

m  en  b,D  ad(r»-l)  aon 
1  P 


=  1  ,  J  +  r-r —  [E  .  -E  ] 

e(n-l)  en  b,D  ad(n-l)  aon 

1  p 


(130) 


Introducing  equations  (121)  through  (124)  into  equation 
(130)  clearing  of  fractions  and  reducing  to  the  simplest  terms, 
the  following  is  obtained: 
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Q  * 


kq  /cf  f  o 

.  Hf.  +  —  [  — 1^  -  f.  +  4  f  j } 

ll  13  1  b.D  P  1  4 

1  p 

/cJc  f 

1+f  (j-l)+  {  ~  [l+f  (j-1)  ]-f  [icf  (j-l)+l]J 

1  p  p  J 


Then  from  equations  (125)  to  (131)  we  find: 


bl(D+V{fl  +  b.D 


[  — 2__i.  _  £  -*■  4f  2 1 } 
l  ^  <*r4  j; 


1  P 


k  ® 
e 


f2  +  b“T  [  'T-  '  K  fs  1 

1  p  M 


(132) 


Since  the  pore  size  of  powder  is  difficult  to  measure 
but  the  particle  size  and  porosity  can  be  measured,  it  is 
desirable  to  substitute  for  Dp  a  function  of  the  porosity,  P, 
and  the  particle  size,  D: 


D 

P  =  -J? 

D+D  (133) 

P 


or: 


D 

P 


PD 

(1-P) 


(134) 


and: 


D  4  0 

P 


D 

(1-P) 


(135) 


Then  the  conductivity  becomes: 


(131) 
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(136) 


D  k  o  (b . 0 PD f  +k  (1-P)  Uf  f-/3f  +4/3 f  2]  j 
1  1  q  6  114 

( 1— P)  fi  ( PD  kaf0  +  2n2  k  (1-P)  [  f  _f  _-/3f  _  2  ]  i 
2  g  2  6  5 

Note  that  the  tc  in.  the  second  part  of  the  denominator 

has  been  removed  and  2b/3/k.a  has  been  introduced  in  its  place, 

(where  k  is  the  thermal  conductivity  of  the  solid)  and  that 

the  equation  has  been  cleared  of  fractions.  At  t?  s  point, 

co  simplify  the  above  expression  to  one  that  can  be  used,  it 

is  necessary  to  make  a  tedious  substitution  for  f  through 

f  from  equations  (85)  through  (91).  This  will  be  done  only 

for  the  case  of  a  vacuum;  i.e. ,  where  k  is  zero. 

y 

The  result  is: 

2/3  D  b^  (1-po)  (l+/0  £>inh  oD 

(1-P)  {2(l-p  )  (coshaD-1) +[  2/3(po  +  Pi)  (l+/c) +<caD(l-pi)  ]sinh  oD] 

(137) 

This  equation  illustrates  the  usefulness  of  the 
general  theory  in  that  a  relatively  simple  equation  can  be 
derived  for  a  very  complex  problem.  This  equation  has  proper 
end  points,  partially  due  to  the  introduction  of  the  surface 
reflections. 


3,2  Validity  of  the  Hamaker  Approximation 

In  development  of  the  Hamaker  equations,  an  approxima¬ 
tion  is  introduced  in  the  conductive  part  of  the  combined  heat 
transfer  by  conduction  and  radiation.  Measurements  of  the 
effective  thermal  conductivity  of  polycrystalline  alumina  at 
temperatures  of  about  2000°K  indicate  that  the  conductive  part 
is  still  a  major  contributor  (^)  u  This  results  from  the 
relatively  high  conductivity  and  large  scattering  in  dense 
polycrystalline  ceramics.  Since  the  utility  of  the  specific 
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solutions  of  the  Hamaker  equation  derived  thus  far  depends 
critically  on  the  validity  of  our  underlying  assumptions,  it 
is  essential  to  have  critical  understanding  of  the  nature 
of  these  limitations. 


3.2.1  The  Hamaker  Approximation 

The  equations  describing  the  transmission  of  heat 
by  conduction  and  radiative  transfer  in  a  homogeneous 
medium  are  (for  the  one  dimensional  case  and  the  steady 
state) : 


k  _  2an2  c*  T4  +  a(I+J)  =  0  (138) 

dx2 

*  an2  o'  T4  -  (a+s) I  +  s  J  (139) 

“  *  (a+s)  J  -  si  -  an2  a'  T4  (140) 


where  I  and  J  are  the  radiant  fluxes  in  the  +x  and  -x  direc¬ 
tions,  T  is  the  absolute  temperature,  a'  is  the  Boltzmann 
constant,  and  k,  a,  s,  n  are  material  constants  already  defined. 
These  equations  are  identical  to  (51)  (55)  and  (56) • 

Hamaker 's  approximation  consists  of  replacing  the 
first  term  in  equation  (138),  k  ^_J£,  by  k  ■-  d2  (g4)  #  in  which 

dx2  To 4  dx2 

To  is  some  average  temperature  in  the  material.  This  replace¬ 
ment  makes  the  set  of  equations  (138) ,  (139) ,  (140)  linear  and 
homogeneous  in  the  quantities  I,  J,  and  T^  so  that  a  solution 
in  closed  form  can  be  obtained  readily.  However,  consider 
the  identity: 
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d2T 

dx^ 


3k 

T 


(f) 

dx 


4T- 


2  4 


dx' 


(141) 


Hamaker ' s  approximation  is  equivalent  to  neglecting  the 
first  term  on  the  right  hand  side:  it  should  be  valid 
therefore  only  when  the  inequality,  (5) , 


<f> 

dx 


« 


T 

3 


d^T 

dx^~ 


(142) 


2  ,  2 

is  satisfied  or  when  the  entire  conduction  term,  k  d  T/dx 
negligible  compared  to  the  other  terms  in  equation  (138). 

We  note  that  the  inequality  (142)  cannot  be  satisfied  with  a 
constant  temperature  gradient  (i.e.  when  conduction  dominates 
radiative  transfer)  and,  in  fact,  is  not  satisfied  for  any^ 
temperature  distribution  of  the  form  T  »  axn  unless  n  «  < 


3.2.2  Numerical  Integration  of  Exact  Equations  at 
Low  Temperatures 


For  computation  it  is  convenient  to  reduce  the  set  of 
equations  (138),  0  39),  (140)  to  a  single  second  order 
equation.  Equations  (139)  and  (140)  are  added  and  the  sum 
integrated  once  to  give: 


I  +  J  -  - (a+2s)  (Hx+kT)  +  A 


(14  3) 


in  which  A  i.s  a  constant  of  integration  and  H  is  the  net 
heat  flux  (also  constant)  defined  by: 


H  =  I  -  J  -  k 


dT 

dx 


(144) 
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Substituting  in  equation  \138) ,  we  obtain  finally: 


K  ”  2an2a ' T4  -  a(a+2s) (Hx+kT)  +  aA  =  0  (145) 

This  equation  has  been  given  essentially  by  Kellett;  we  have 
introduced  here  the  effect  of  radiation  scattering  and  multi¬ 
plied  the  Boltzmann  constant  by  n^.  For  comparison  with 
various  approximations  we  have  calculated  an  exact  solution 
to  equation  (145)  by  numerical  integration  for  the  initial 
conditions  (imposed  at  x=0)  listed  below: 

a.  Heat  flow  in  -x  direction  in  a  homogeneous 
medium  containing  the  plane  x=0 

b.  Net  flux,  H  -  -3.03  cal.  cm  2  sec  ^ 

c.  Radiant  flux  in  +x  direction,  I(x)  =  0  at  x  -  0 

d.  Radiant  flux  in  -x  direction,  J(x)  -  -0.03  at 
x=0  (cal  cm"2  sec*”!) 

e.  Temperature,  T(x)  =  0.1  at  x~0  (in  units  of 
1000°K) . 

dT 

f.  Temperature  gradient,  ~  -  1.0  at  x=0  (units  of 

1000 °K  cro~i)  dx 

The  matei*ial  constants  assumed  are: 

g.  Absorption  coefficient,  a  =  10  cm  1 

h„  Scattering  coefficient,  s  -  10  cm  1 

-1  -i 

i.  Material  conductivity,  k  -  3  cal  cm  sec  per 
1000 °K 

j.  Refractive  index  and  radiation  constant, 
n2i'  =3  cal  cm  2  sec  ^  per  (1000°K)^ 
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For  the  particular  values  assumed,  equation  (145) 
becomes: 


=  20  T4  +  300  T  -  303  x  -  30.1  (146) 

dx^ 


A  straightforward  numerical  integration  gives  the  solutions 
plotted  in  Figure  1.  It  should  be  noted  that  the  initial 
temperature  was  chosen  to  be  very  low  and  the  temperature 
gradient  is  large  so  the  solutions  obtained  thus  far  cover 
a  range  in  which  substantial  differences  between  the  Hamaker 
and  exact  solutions  are  to  be  expected. 

3.2.3  Numerical  Integration  of  Exact  Equations  at 
High  Temperatures 


We  have  also  attempted  the  numerical  integration  of 
the  transfer  equation  (145)  in  the  temperature  range  of 
1000°K  and  above.  In  principle  a  numerical  solution  can 
be  obtained  co  arbitrarily  high  accuracy  by  taking  sufficiently 
small  intervals  for  the  integration?  in  practice  a  very  large 
number  of  intervals  is  needed  to  get  reasonable  accuracy  in 
tie  high  temperature  region.  As  the  solution  proceeds  to 
higher  temperatures  the  second  term,  2  an2c'T"*,  dominates; 

d2T 

to  satisfy  the  equation  the  first  term,  k  ^  2  »  must  become 
large  and  positive  with  the  result  that  the  solution  T(x) 
rapidly  goes  to  infinity. 

A  complete  analytical  solution  for  equation  (145)  has 
not  been  found  but  we  have  developed  a  form  which  has  the 
correct  asymptotic  behavior  and  is  well  adapted  to  numerical 
calculation  in  the  high  temper at  ire  limit.  We  can  show 
that  if  T(x)  is  a  solution  of  equation  (145) ,  then  the  inverse 
function  x(T)  must  satisfy  equation  (147): 


Heat  flux  ( ca!  sec~!cm-2) 


x  (T) 


(05-To5)  +  a(a+2s)(G2-T02) 


2  , 
an  o  1 

k 


2aA 
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(e-T0)  + 


dT 

(--)  + 

QXo 


2a (a+2s) 
k 


a  1 

H  /  x(^)dej"  2 
To 


(147) 


The  two  variables  of  integration,  9  and  have  the 
dimensions  of  temperature:  0  ranges  from  T0  to  T  and  £ 
ranges  from  T0  t~*  0. 

After  differentiating  once ,  equation  (147)  becomes: 


^  „  2  , 

dx  _  r4_  an  a 1 

dT  ~  *  5  k 


(T5-T05)+a(a+2s)  (T2-T0 2) 


2aA 

k 


(T-To )  + 


/dT.  2  2a(a+2s) 

dXo  k 


T  -  * 

H  /  x  (£ )  d  £;  }  2 

To 


(148) 


Although  formidable  at  first  sight,  these  equations 
are  not  difficult  to  apply.  A  procedure  for  calculating 
the  solution  x  (T)  follows: 


a. 


b. 


c. 


Assume  the  materials  constants,  the  heat  fluxes  H 
and  A,  and  the  initial  conditions  T0  and(^)  at  xo=0 

Estimate  the  coordinate  x.  corresponding  to  a  nearby 
temperature  T  and  evaluate  the  integral 

T1 

/  x  a)  de. 

To 

(Jv 

Calculate  ~  at  T  =  T^  from  equation  (148). 


T 

1  dx 

d.  Integrate  to  get  x  =  /  (~ )  dT  +  x0. 

1  To 

e.  Use  the  calculated  value  of  x^  to  replace  the  esti¬ 
mated  value  in  (b)  and  iterate  until  a  consistent 

value  for  x,  is  obtained. 

1 

f.  Continue  stepwise  to  calculate  x2  (T2) ,  x^  (T^)  and 

so  forth.  As  the  solution  proceeds  to  high  tempera¬ 
ture  the  step  (b)  and  the  iteration  becomes  unnecessary 

T 

*  2^ 

because  the  term  involving  /  x(£,)d£  becomes 

To 

negligible  compared  to  the  first  term  on  the  right 
hand  side  of  equation  (148). 


As  an  illustration  of  the  method  we  have  calculated 
a  solution  with  the  same  material  constants  and  initial  values 

dT 

for  H,  I,  J,  and  (“)  used  in  the  last  section  except  that  the 
initial  temperature  °  is  taken  to  be  1000°K  (To=1.0) .  Equation 
(148)  then  becomes: 


dx 

dT 


(8(T5-1)  +300(T2-1)  -60. 2(T-l)+l-606  /T  x  d 

1 


1 

r  2 

(149) 


The  result  is  plotted  in  Figure  2. 

To  show  that  equations  (145)  and  (147)  are  equivalent, 

differentiate  equation  (147)  with  respect  to  T  to 
obtain  equation  (148)  which  can  be  written  as: 

(dT}  =  {  5  (t5-t°5}  *  a(a+2s)  (T2-To2)  - 

if  lT-T°)  +  (S>2  +  H  s'.  (O  *)  *  (150) 

To 
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